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Abstract 

The strength distribution of the isoscalar giant dipole resonance (ISGDR) in 58 Ni has been 
obtained over the energy range 10.5-49.5 MeV via extreme forward angle scattering (including 0°) 
of 386 MeV a particles. We observe a "bi-modal" El strength distribution for the first time in 
an A < 90 nucleus. The observed ISGDR strength distribution is in good agreement with the 
predictions of a recent RPA calculation. 



PACS numbers: 24.30.Cz, 21.65.+f, 25.55.Ci, 27.40+z 
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The compressional-mode giant resonances in the atomic nuclei — the isoscalar giant 
monopole resonance (ISGMR) and the isoscalar giant dipole resonance (ISGDR) — provide 
a direct method to obtain the incompressibility of the nucleus and of nuclear matter (K nm ) 
[lj]. Although ISGMR has been investigated extensively for a large number of nuclei in the 
past, the exotic ISGDR has been identified only in a few nuclei and the location of ISGDR is 
not systematically established over the wide mass region. One major concern with ISGDR 
data had been that the nuclear incompressibility extracted from the centroid of the ISGDR 
strength distribution was significantly different from that obtained from the ISGMR data. 
In recent work, this ambiguity has been resolved for 208 Pb by a more precise, background- 
free measurement of ISGDR strength distribution, and the value of K nm obtained from the 
ISGMR data is now consistent with that from the ISGDR data for 208 Pb Q- 

The experimentally-observed ISGDR strength distribution in all A > 90 nuclei has a 
"bi-modal" structure ^, U, 0, 0], in agreement with predictions of recent theoretical work 
. Of these, only the high-energy (HE) component depends on K nm and, hence, 
is of interest from the point of view of determining an experimental value for this important 
parameter. The low-energy (LE) component, which is quite small in comparison with the HE 
component, is located much higher in excitation energ y th an the expected Ifkv component of 
the ISGDR, previously identified by Poelhekken et al. [11]. As well, it is lower in energy than 
the isovector giant dipole resonance (IVGDR) which can be excited in inelastic a scattering 
via Coulomb excitation; in the event, the full expected IVGDR strength is subtracted out 
in the analysis of all aforementioned data. The exact nature of this component is not fully 
understood yet, although sugges tions have been made that it might represent the "toroidal" 
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131 ] provide a review of the recent experimental and theoretical 



or "vortex" modes; Refs. 
results on ISGDR. 

For 58 Ni, there has been only one recent measurement, wherein a concentrated ISGMR 
and isoscalar giant quadrupole resonance (ISGQR) strength distribution has been observed, 
but the ISGDR strength is reported to be spread more or less uniformly over E x = 12 to 35 
MeV j^]. This observation leaves a few open questions: Is the ISGDR strength fragmented 
in light nuclei such as 58 Ni? Are we missing the resonance strength distribution because 
of experimental limitations? In an attempt to answer these questions, we have carried out 
measurements on excitation of isoscalar giant resonances in 58 Ni. In this Letter, we report 
our results on the ISGDR strength distribution in 58 Ni. We find that the ISGDR in this 
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nucleus has a "bi-modal" structure as well, similar to that in the medium- and heavy-mass 
nuclei, and that the experimentally observed ISGDR strength is in good agreement with 
predictions of a recent RPA calculation. 

The 58 Ni(a, a') experiment at E a = 386 MeV was performed at the ring-cyclotron facility 
of Research Center for Nuclear Physics (RCNP), Osaka University. Details of the experi- 
mental measurements and data analysis procedures have been provided in Refs. j^Q]; only 
the salient points are elaborated upon below, a-particles, inelastically scattered off a 5.8 
mg/ cm 2 -thick 58 Ni target, were momentum analyzed in the spectrometer, Grand Raiden 
and detected in the focal-plane detector system comprised of two multi-wire drift- 
chambers and two scintillators, providing particle identification as well as the trajectories of 
the scattered particles. The scattering angle at the target and the momentum of a scattered 
particle were determined by the ray-tracing method. The 58 Ni(a, a') spectra were measured 
in the angular range of 0° to 8.5° for two excitation-energy-bite settings of the spectrometer 
(E x = 5.0-35.0 MeV and E x = 22.0-52.0 MeV). The primary beam was stopped at one of 
four different Faraday cups, depending on the scattering angle and the excitation energy 
bite of the spectrometer. The vertical position spectrum obtained in the double-focused 
mode of the spectrometer was exploited to eliminate the instrumental background due to 
Coulomb scattering of the beam at the target and subsequent rescattering by the edges of 
the entrance slit, the yoke, and walls of the spectrometer 0, Fig. 1 shows an excita- 
tion energy spectrum for the 58 Ni(a, a') reaction at 9 avg , = 0.69° after subtraction of the 
instrumental background. A prominent "bump" corresponding to (ISGMR + ISGQR) in 
58 Ni is observed at E x = 10-25 MeV and another bump [ISGDR + the high-energy oc- 
tupole resonance (HEOR)] is visible as a shoulder at E x ~ 33 MeV. There is an underlying 
continuum in the high excitation-energy region in the spectrum. Since there is no sound 
theoretical basis to estimate and subtract the physical continuum from the excitation en- 
ergy spectrum, it is reasonable to assume that the continuum background remaining after 
elimination of the instrumental background is the contribution from the higher multipoles 
and the three-body channels resulting, for example, from knock-out reaction. In the present 
work, a multipole-decomposition (MD) analysis has been performed to extract giant reso- 
nance strengths, by taking into account the transferred angular momentum up to AL = 
7. The cross-section data were binned in 1-MeV energy intervals to reduce the statistical 
fluctuations. For each excitation-energy bin from 10.5 MeV to 49.5 MeV, the experimental 
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angular distribution <J exp {9 c . m ., E x ) has been fitted by means of the least-square method with 
the linear combination of calculated distributions cr cal (6 c .m., E x ) defined by: 



L=7 



a exp (8 c . m .,E x ) = £ a L (E x ) x af(9 c . m ., E x ), (1) 



L=0 



where cr|^(0 cm , E x ) is the calculated distorted-wave Born approximation (DWBA) cross 
section corresponding to 100% energy-weighted sum rule (EWSR) for the L th multipole. 

The DWBA calculations were performed following the method of Satchler and Khoa 
[if] ] using density- dependent single folding, with a Gaussian a-nucleon potential (range 
t = 1.88 fm) for the real part, and a Woods-Saxon imaginary term; the calculations were 



carried out with the computer code PTOLEMY 



171 ] . Input parameters for PTOLEMY were 



modified ^| to take into account the correct relativistic kinematics. The shape of the real 
part of the potential and the form factors for PTOLEMY were obtained using the codes 
SDOLFIN and DOLFIN Q]. We used the transition densities and sum rules for various 
multipolarities described in Refs. HQ. The radial moments for 58 Ni were obtained by 
numerical integration of the Fermi mass distribution with c = 4.08 fm and a = 0.515 fm 
[20I 1 . The folding- model parameters with the computer code PTOLEMY were obtained from 
analysis of 58 Ni+ a elastic- and J n = 2 + inelastic-scattering data at E a = 386 MeV taken 
in a separate experiment. The folding model parameter extracted for the real part of the 
potential is V = 37.02 MeV, and the parameters for the Woods-Saxon type imaginary part 
were: W = 36.86 MeV, r/(reduced radius) = 0.95 fm, and a/(diffuseness) = 0.67 fm. Using 
these parameters, the DWBA calculation for the first J 71 " = 2 + state in 58 Ni was carried out 



with 

e 2 b 2 



TOLEMY using a collective form factor with the previously-known B(E2) = 0.070 
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Fig. 2 compares the results of the calculations and the experimental data; 
the calculations reproduce elastic scattering cross sections as well as the inelastic scattering 
differential cross section for the 2 + state very well. 

The contribution of the IVGDR excitation to the measured cross sections was subtracted 
prior to multipole decomposition. The cross section for IVGDR excitation was calculated 
using the strength distribution obtained from photonuclear work (2^ in conjunction with 
DWBA calculations on the basis of the Goldhaber- Teller model. The fits of the angular 
distributions for two energy bins near the peaks of the ISGMR and ISGDR are shown in 
Fig. 3; the L — 0, 1, 2 and 3 contributions to the differential cross section are also shown. 
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The ISGMR has a maximum at 9 c , m , = 0° and its contribution is dominant in comparison to 
the other multipoles at E x = 18.5 MeV. Similarly, the ISGDR has a maximum at 9 c _ m , ~ 2° 
and its contribution is dominant in comparison to the other multipoles at 29.5 MeV. The 
fitted parameters ai{E x ) so obtained are fractions of the EWSR's, which can be related to 
the strengths Sl(E x ) as follows: 

S (E X ) = 2h2A< J 2> a (E x ), (2) 
mE x 

S l (E x ) = -^4_(n < r 4 > < r 2 > 2 -lOe < r 2 >) ai (E x ), (3) 

S> 2 (E X ) = 7^rH2L + l) 2 < r 2L ~ 2 > a h (E x ), (4) 



where m, A, and < r N > are the nucleon mass, the mass number, the N t/l moment of the 

■2 



ground-state density, respectively, and e=(4/E 2 +5/E )fr /3mA. E and E 2 are the centroid 
energies of the GMR and GQR, respectively. The strength distributions extracted from 
these fits for L = (ISGMR), L = 1 (ISGDR), and L = 2 (ISGQR) in 58 Ni are shown 
in Fig. 4. In order to examine the reliability of the strength distributions obtained from 
the fits in the MD analysis, we varied the L max value from L = 6 to L = 8. However, 
the extracted strength distributions for L— 0-3 did not change in any significant way. In 
addition, a completely independent data analysis, using a different folding-model potential, 
led to essentially the same results for the various strength distributions. 

The centroid energy of ISGMR, shown in Fig. 4(a), was determined to be E x = 19. 
MeV between E x = 10.5 and 32.5 MeV. A total of 92t 4 % of the E0 EWSR was identified 
in the above excitation-energy range. [The errors quoted in all EWSR values here are only 
statistical; in addition, there may be a 15-20% systematic error in the EWSR fractions 
because of the uncertainties associated with the DWBA calculations used in the MDA anal- 
ysis.] This result is similar to that reported earlier where the fraction of 74t 2 2% for the 
E0 EWSR value is observed between E x = 12.0 to 31.1 MeV with a centroid of 20.30to;?4 
MeV. 

The strength distribution of ISGDR is shown in Fig. 4(b). We observe a "bi-modal" 
strength distribution between E x = 10.5 and 42.5 MeV. A low-energy (LE) component at 
~16 MeV appears as a shoulder at the low-energy side of the extracted ISGDR strength. The 
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excitation energy of this component is much higher than that expected for the previously- 



mentioned lhu component of the ISGDR 11], and lower than that of the IVGDR (E x = 18 
MeV This is the first observation of concentrated "bi-modal" isoscalar El resonance 

in an A < 90 nucleus. In the previous measurement on ISGDR strength distribution in 
58 Ni [1J], it was reported that 41% El EWSR was spread more or less uniformly over E x 
= 12 to 35 MeV. In the present case, not only is the ISGDR strength concentrated and has 
a "bi-modal" distribution similar to that observed in the A > 90 nuclei, it would appear 
that nearly all of the expected ISGDR strength has been observed: A total of 98^5% of 
the El EWSR was identified between E x = 20.5 and 40.5 MeV. The centroid energy of 
the ISGDR for the same excitation energy range is determined to be 30.8±^ MeV. The 
difference between the present result and the result of Ref. Q| might be attributable to 
the fact that in Ref. |l4| the multipole decomposition is carried out after subtracting a 
"background" from the excitation energy spectrum, whereas, as pointed out earlier, no such 
subtraction is required in the present analysis since our spectra have been rendered free of 
all instrumental background. 

The ISGQR strength distribution is shown in Fig. 4(c). In this case, 73±|% of the E2 
EWSR value has been observed in the excitation energy range between E x = 10.5-21.5 MeV. 
The centroid energy of ISGQR is determined to be 16.3to g. This result is consistent with 
the result of Ref. where E2 strength corresponding to 115±18% of the E2 EWSR 

was found between E x = 10.5-20.5 MeV with a centroid of 16.1±0.3 MeV and rms width 
of 2.4±0.2 MeV. The ISGQR strength shows a near constant value beyond E x = 20 MeV. 
At present, the reasons behind this extra strength are not fully understood. Similarly, en- 
lanced El strengths at high excitation energies have been noted previously in other nuclei 
^,0, and have been attributed to contributions to the continuum from three-body channels, 
such as knock-out reactions. These processes are only implicitly included in the MD anal- 
ysis as background and may lead to spurious contribution to the extracted giant resonance 
strengths at higher energies where the associated cross sections are very small. This con- 
jecture is supported by recent charged-particle decay measurements on ISGDR wherein no 



such spurious strength at high excitation energies is observed 



12, 



24]. Incidentally, a similar 
increase at higher excitation energies has also been reported recently in E0 strength in 12 C 
when a multipole decomposition was carried out without subtracting the continuum 
from the excitation-energy spectra. 
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To get a quantitative understanding of ISGDR strength distribution, the experimental 
strength distribution has been compared with the predictions of quasi-particle random phase 
approximation (QRPA) as shown in Fig. 5. The agreement between the experimental 
and theoretical ISGDR strength distributions is rather good, except at the very highest 
excitation energies where, as mentioned previously, the experimentally extracted ISGDR 
strength is compromised by the limitations of the MD analysis procedure. This is quite 
remarkable since, in general, the details of the theoretical strength distributions do not 
quite match the experiment even when the centroid energies are in agreement with the 
experimental data. 

In summary, we have performed 58 Ni(a, a') measurements at E a = 386 MeV to study the 
excitation of ISGDR. The ISGDR strength distribution has been obtained up to E x = 49.5 
MeV by multipole-decomposition analysis. A two-component ISGDR strength distribution 
has been observed for the first time in 58 Ni and, indeed, in any A < 90 nucleus. The centroid 
energy of the high-energy component of the ISGDR (E x = 30.81^1 MeV) is consistent with 
the global systematics and the strength distribution for the HE component of the ISGDR is 
in good agreement, qualitatively and quantitatively, with predictions of QRPA calculations 
with a nuclear matter incompressibility value of K nm = 217 MeV 

This work has been supported in part by the US- Japan Cooperative Science Program of 
JSPS and the US National Science Foundation (grants number INT03-42942 and PHY04- 
57120). 
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FIG. 1: Excitation energy spectrum for the 58 Ni(a,a') reaction at E a = 386 MeV. Inelastically- 
scattered a particles were measured with the magnetic spectrometer at 9 = 0° with two different 
settings of the magnetic field to cover the excitation-energy ranges of £^=5.0-35.0 MeV and E x = 
22.0-52.0 MeV. 
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FIG. 2: Cross sections for (a) the a + 58 Ni elastic scattering and (b) the 58 Ni(a, a') 58 Ni(2 + ) 
reaction at E Q = 386 MeV. The solid lines are the results of the folding model calculations. 
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FIG. 3: Angular distributions of selected 1 MeV bins for the 58 Ni(a,a') reaction at 386 MeV. 
(a) Results for E x = 18.5 MeV. The open circles are the experimental data. The lines show 
contributions from L = (dotted line), L = 1 (dashed line), L = 2 (dot-dashed line), L = 3 
(double dot-dashed line) and other higher-multipole components including IVGDR (double dash- 
dotted line), respectively, (b) Same as part (a), except for E x = 29.5 MeV. 
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FIG. 4: Strength distributions for the L = 0, 1, and 2 in 58 Ni. (a) ISGMR, (b) ISGDR, (c) 
ISGQR. The errors shown for each excitation-energy bin were estimated by changing the strength 
parameter for one component in order to satisfy A% 2 increase by 20% while fitting with the other 
parameters remaining free. The low-energy component of the ISGDR is indicated by an arrow. 
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